Thermodynamic and thermoelectric properties of (Ga,Mn)As and related compounds 
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Various experimental results providing information on thermodynamic density of states in 
(Ga,Mn)As are analyzed theoretically assuming that holes occupy GaAs-like valence bands. Al- 
lowing for Gaussian fluctuations of magnetization, the employed model describes correctly a critical 
behavior of magnetic specific heat found experimentally in (Ga,Mn)As near the Curie temperature 
Tc [S. Yuldashev et al., Appl. Phys. Express 3, 073005 (2010)]. The magnitudes of room temper- 
ature thermoelectric power, as measured for GaAs:Be and (Ga,Mn)As [M. A. Mayer et al., Phys. 
Rev. B 81, 045205 (2010)], are consistent with the model for the expected energy dependencies of 
the hole mobility. The same approach describes also temperature variations of conductance specific 
to the Anderson-Mott localization, found for various dimensionality (Ga,Mn)As nanostructures at 
subkelvin temperatures [D. Neumaier et al, Phys. Rev. Lett. 103, 087203 (2009)]. We conclude 
that the examined phenomena do not provide evidence for an enhancement of density of states by 
the presence of an impurity band at the Fermi energy in ferromagnetic (Ga,Mn)As. Furthermore, we 
provide for (Ga,Mn)As expected values of both electronic specific heat at low temperatures T< Tc 
and magnetization as a function of the magnetic field at Tc- 
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I. INTRODUCTION 

Extensive studies of hole-controlled ferromagnetic 
semiconductors, such as (Ga,Mn)As and p-(Cd,Mn)Tc, 
have allowed to demonstrate a variety of functionali- 
ties specific to these systems Jul At the same time, how- 
ever, results of various experiments have faced us with 
a number of challenges suggesting that the understand- 
ing of these materials is by far not satisfactory. For 
instance, recent studies of (Ga,Mn)As using two kinds 
of tunneling spectroscopy have lead to two entirely con- 
tradictory pictures of electronic states in the vicinity of 
the Fermi energy Ep. According to scanning tunneling 
microscopy^ 3 - and in agreement with the Altshuler and 
Aronov 4 and Finkelstein- description of the Anderson- 
Mott localization for the spin polarized band carriers, 6 
the one-electron density of states (DOS) attains a mini- 
mum at E-p and its depression extends over an energy 
range of the order of the momentum relaxation rate, 
Tl/t ~ 100 meV in (Ga,Mn)As/ In contrast to one- 
electron DOS probed in tunneling experiments, the DOS 
for charge excitations as well as thermodynamic DOS, 
Py = dp/dEp, are only weakly renormalized by carrier 
correlation and disorder in doped semiconductors on the 
metal side of the Anderson-Mott transition, where the 
ratio of the inter-particle distance to an effective Bohr 
radius is relatively small r s — (47rp/3)~ 1 / 3 /a| 3 < 2.4. 

A rather different picture emerges from resonant tun- 
neling spectroscopy of (Ga,Mn)As quantum wells4^ Ac- 
cording to the interpretation of the accumulated findings, 
the Fermi level resides in a narrow impurity band lo- 
cated ~ 50 meV above the edge of hole subbands. These 
subbands are well resolved by resonant tunneling spec- 
troscopy and, thus, virtually unaffected by disorder^ 

The impurity band model implies an enormous ther- 



modynamic DOS of carriers at the Fermi level. An un- 
usually large effective mass emerges also within some 



interpretations of the optical^ and transport 11 proper- 
ties of (Ga,Mn)As. Historically, large magnitudes of ef- 
fective masses in, for example, heavy-fermion systems^ 
and Kondo alloy a 13 ' 14 were discovered by studies of ther- 
modynamic and thermoelectric effects. Therefore, it is 
of particular importance to describe theoretically those 
properties of (Ga,Mn)As, which provide information on 
thermodynamic DOS. Since this DOS, in contrast to one- 
electron DOS probed in tunneling experiments, does not 
exhibit a Coulomb anomaly at EpJ> it should show a 
significant enhancement if relevant electronic states are 
indeed characterized by a large effective mass. 

In this paper, various experimental results providing 
information on thermodynamic DOS in (Ga,Mn)As are 
analyzed theoretically assuming that holes occupy GaAs- 
like valence bands, described here by the six band k ■ p 
model with effects of the p-d exchange interaction in- 
cluded within the molecular field approximation! 15 i 16 We 
show that the employed model, taking into account Gaus- 
sian fluctuations of magnetization, describes correctly 
a critical behavior of magnetic specific heat found ex- 
perimentally in (Ga,Mn)As around the Curie temper- 
ature TqJ^ We also provide the expected field depen- 
dence of magnetization at Tq within the mean-field ap- 
proximation. Furthermore, we show that the magnitudes 
of room temperature thermoelectric power, as measured 
for GaAs:Be and (Ga,Mn)As^ are consistent with our 
model for the expected energy dependencies of the hole 
mobility. We then turn to low temperatures T <C Tc 
and present computed values of the Sommerfeld coef- 
ficient 7, describing electronic specific heat, and mag- 
netization as a function of the magnetic field at Tq. 
We also employ the same approach to discuss temper- 
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ature variations of conductance specific to the Anderson- 
Mott localization, determined experimentally for various 
dimensionality (Ga,Mn)As nanostructures at subkelvin 
temperatures. 19 We conclude that the examined phenom- 
ena do not provide any evidence for an enhancement 
of thermodynamic DOS by the presence of an impurity 
band at Ep in ferromagnetic (Ga,Mn)As. 



II. CRITICAL BEHAVIOR OF SPECIFIC HEAT 

Recently, a critical behavior of specific heat in 
Ga 1 _ 2 ,Mn :I .As was resolved experimentally for two sam- 
ples with the Mn concentration x = 1.6 and 2.6%, which 
showed insulating and metallic behavior, respectively^- 7 - 

In order to describe these findings we employ the 
Ginzburg-Landau approach taking into account critical 
fluctuations in the Gaussian approximation^ We spec- 
ify by (Ti(x), i = 1,2, ... ,n, a local magnitude of an n- 
dimensional order parameter (e. g., a magnitude of local 
spin density) around x in a d dimensional block. The 
block Hamiltonian (the free energy functional) H[a(x)] 
at temperature T and in the magnetic field h can then be 
written as an expansion in powers of a{x) and Wa{x)^ 



H[a{x)]/k B T = Jd d x [a +a 2 o 2 



(1) 



+ a 4 (a 2 ) 2 +c(Va) 2 -h-a], 



where 



a 2 = S ^{p l {x)f , 

i=l 

d n / a \ 2 
a — 1 i—l v 7 



(2) 
(3) 



The temperature dependent coefficients (k = 0,2,4) 
and c describe the free energy cost associated with a 
change in the magnitude and in the local direction of 
the order parameter, respectively. 

Within the Gaussian approximation and taking into 
account only terms that are singular on approaching the 
ordering temperature Tq , the specific heat assumes a crit- 
ical behavior given bj*2£ C = C ± i d / 2-2 , where C + = nCa 
and C~ = 2 rf / 2 C apply to the T > T c and T < T c cases, 
respectively; t = \T — Tq\/Tc', and 



Co 



(27r)-V 2 T c /c)2 / d d k'(l + k' 2 ) 



(4) 



where a' 2 defined by a 2 = a' 2 (T — Tc) shows no singular 
dependence on temperature around To 

In order to determine the magnitudes of a' 2 and c 
for (Ga,Mn)As we note that, within the p-d Zcner 
model ; 15 ' 21 ' 22 the free energy functional consists of two 
contributions. The first describes the free energy of lo- 
calized spins of magnitude S in the absence of carriers, 
whereas the second is the free energy of the Fermi liq- 
uid of holes in the valence band in the molecular field 



of a prescribed configuration of the localized spins (for 
such a description to be possible it is required that the 
dynamics of the localized spins are much slower than 
that of the carriers). The form of particular contribu- 
tions to the Ginzburg-Landau free energy functional was 
determine d 15 ' 22 ! 23 employing magnetization of localized 
spins M(x) = a(x)Msat/S as an order parameter, where 
the saturation magnetization is related to the magnetic 
moment g^S and the effective concentration of local- 
ized spins Nqx c s in a standard way Mg at = g^BSN x e s, 
where the g-factor g m 2.0. Neglecting interactions be- 
tween localized spins in the absence of carriers, and as- 
suming that the carrier liquid is strongly degenerate ; 15 ' 22 
we obtain 



a' 2 T c = 



3x eS N Q 



2 " 2S{S+1) 
and, in terms of the magnetic stiffness ^4 ; 23 ' 24 



A{M)M l 
m™o S 2 M 2 k B T c 



Sat 



(N Q p XcS ) 2 B 
k B T c 2 ' 



(5) 



(6) 



where B is a property of the electronic subsystem, Ao is 
the cation concentration and /3 is the p-d exchange con- 
stant. 

In Ref. [13, experimental data were presented for the 
specific heat of two (Ga,Mn)As films. According to x- 
ray diffraction measurements and x-ray microanalysis, 
the Mn concentration x is 1.6% in the sample A and 
2.6% in the sample B, with Tq of 40 and 52 K, respec- 
tively. Due to relatively low x values, these samples 
are expected to be close to the metal-insulator transition 
(MIT) . This expectation is confirmed by temperature de- 
pendencies of resistance, which indicate that samples A 
and B are on the insulating and metallic side of the MIT, 
respectively. 17 According to magnitudes of the Hall resis- 
tance /Oh at room temperature, the hole concentrations 
are 2.7 x 10 19 and 4.5 x 10 19 cm~ 3 for samples A and 
B, respectively. Since, however, pn tends to diverge near 
the MIT— (if measured at temperatures below the impu- 
rity binding energy, ~ 1000 K in GaAs:Mn), it leads to 
underestimated values of the hole concentrations. Fur- 
thermore, owing to critical fluctuations in the local DOS 
near the MIT,— ' 15 ' 26 only a part of the sample volume 
is ferromagnetic! 15 ' 26 This means that the apparent val- 
ues of low-temperature spontaneous magnetization are 
smaller than the expected magnitudes of saturation mag- 
netization Msat = g^BSN Xcs, as observed) 18 ' 27 ' 28 

Assuming that a; c ff < x only due to the presence of in- 
terstitial Mn, we have a; e ff = x — 2x{, where x and x\ are 
the total and interstitial Mn contents. Similarly, the hole 
concentration is p — No(x — 3xi). 29 Using these equations 
and the mean-field theory^i 5 - we determine the values of 
x e s and p, which reproduce the experimental values of 
T c . We find x cS = 1.6% and p = 3.5 x 10 20 cm- 3 for 
the sample A, whereas for the sample B the values are 
2.0% and 3.8 x 10 20 cm -3 , respectively. As expected, the 
magnitudes of the hole concentrations obtained in this 
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FIG. 2. (Color online) Computed values of the proportionality 
coefficient between M(H) 3 and H for Mn magnetization M 
in the magnetic field H at Tc- The calculation has been 
performed within the p-d Zener model for various effective 
Mn contents x in (Ga,Mn)As. 



lated magnitude of the specific heat constitutes an upper 
limit of the expected values. 



FIG. I. (Color online) Theoretical temperature dependence 
of the magnetic specific heat calculated with no adjustable 
parameters for the sample A (upper panel) and the sample B 
(lower panel). The solid and dashed lines correspond to the 
Heisenberg (n = 3) and Ising (n = 1) models, respectively. 
Dots represent experimental data (after Yuldashev et 0,1*2.) . 



way are larger than the ones obtained from the Hall resis- 
tance, as quoted above. Similarly, the measured values of 
spontaneous magnetization 17 are by a factor 4.7 and 4.0 
smaller than Ms a t calculated for the values of x e g for 
samples A and B, respectively. We note that the con- 
tribution of holes' magnetic moment to experimentally 
available magnetization may account for its reduction by 
less than 20%^ 

We use the above sample parameters as well as GaAs 
Luttinger parameters and Nq/3 = —1.2 eV i 15 ' 16 The val- 
ues of B have been determined within the six-band k ■ p- 
modelj 23,24 in which we have neglected . The nu- 
merical value is B — 0.17eV _1 nm _1 for both samples. 
These numbers yield Co = 0.13 J/(mol • K) for the sam- 
ple A and C Q = 0.14 J/(mol • K) for the sample B. The 
corresponding dependencies of the specific-heat on tem- 
perature are shown in Fi g. HI where experimental data 
from Figs. 3 and 4 of Ref. [l7| are depicted with dots. We 
see that the theory describes quite reasonably the experi- 
mental values, particularly if one can assume that, owing 
to magnetic anisotropy, the system is effectively Ising- 
like. Actually, since near the MIT only a part of the 
spins contribute to the ferromagnetic order, the calcu- 



III. CRITICAL BEHAVIOR OF 
MAGNETIZATION 

Within the Gaussian approximation, the temperature 
and field dependencies of magnetization M(T, H) are 
given by the mean-field formula. Since theoretically ex- 
pected values of spontaneous magnetization M(T) were 
already presented we discuss here only M(H) at 
T = Tq. This dependence is determined by the coeffi- 
cient d4 in Eq. [1] according to 

d = (V4a 4 ) 1/3 . (7) 

Within the p-d Zener model = a^g) + 04(c) j where 

9(2S 2 + 2S + 1) 



04(S) 



405 3 (5 + l) 3 



-A^OXeff M 6.2 X lO^NoXeB, (8) 



is the localized-spin contribution and ci4( c ), the carrier 
contribution, is calculated numerically from the expan- 
sion of the carrier free energy in the spin splitting param- 
eter B G = A F 0M/(6g(j, B ), 



F c (B G )=F c (0)-p 2 B 2 G + P4 B 4 G , 



(9) 



where is the Landau parameter describing the 
correlation-induced enhancement of the carrier spin sus- 
ceptibility, A-p — 1.2i 15 ' 32 i 33 , whereas the expansion coef- 
ficients pi are determined by the band structure param- 
eters and the hole concentration. We obtain 
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4fc B Tb(a4(S) + 04(c)) 



(10) 
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FIG. 3. (Color online) Comparison of theoretical (solid line) 
and experimental (symbols) field dependence of magnetiza- 
tion at Tc for the sample A of Yuldashev et alM- The dashed 
line shows theoretical values reduced by a factor of 2.45 to 
match experimental data. 



where 



and 



k B T c = A F p 2 N x eS S{S + l)p 2 /54, (11) 



a 4(c ) = A F (pN Q x cS /6) 4 P 4/k B T. (12) 

The computed magnitudes of the slope M (H) 3 /H at var- 
ious hole concentrations and effective Mn contents are 
presented in Fig. [2] The decrease of the slope with in- 
creasing hole concentration, seen in the plot for low x, 
is due to an increase of To The value for the sample A 
is 0.039 (emu/cm 3 ) 3 /Oe, which corresponds to the solid 
line shown in Fig. [3] against the experimental dataJi As 
seen, the theoretical values of magnetization M(H) re- 
produce satisfactorily the character of the field depen- 
dence observed experimentally. However, the absolute 
magnitudes of computed magnetization are by a factor 
of about 2.45 greater than the experimental ones. We 
attribute this discrepancy to a reduction of the ferro- 
magnetic phase volume by critical fluctuations in the lo- 
cal DOS in the vicinity of the MIT, as discussed in the 
previous section. In fact, this factor is even smaller than 
the ratio of saturated and low-temperature spontaneous 
magnetization, determined to be 4.7 for the sample in 
question, as discussed in the previous section. 



IV. HIGH TEMPERATURE 
THERMOELECTRIC POWER 

In general, thermoelectric power contains diffusion and 
phonon drag contributions, 3 - as well as, in the ferro- 



magnetic case, a magnon drag term. The magnitudes 
of the drag terms scale with phonon and magnon relax- 
ation times, so that they dominate at low temperatures, 
particularly in annealed (Ga,Mn)As samples* 3 ^— where 
a reduced concentration of Mn interstitials may suppress 
relevant scattering. 

Awaiting for a quantitative theory of phonon and 
magnon scattering in (Ga,Mn)As, we limit our consid- 
erations to high temperatures, where the diffusion term 
is expected to dominate. For spherical bands with arbi- 
trary dispersion, the diffusion thermopower for carriers 
with charge +e is given by 3 -^ 



S 



where /q in 



dE 



dfo 
dE 



A(E)k 3 (E) 



(13) 



(14) 



is the Fermi-Dirac distribution function and k 3 (E) de- 
scribes the spherical band. The mobility fi depends on 
the hole energy E with respect to the top of the corre- 
sponding valence band subband as 



H(E) = voE r , 



(15) 



where the exponent r depends on the mechanism which 
limits the carrier mobility. We have r w —1/2 for scat- 
tering by acoustic phonons and r rs 3/2 for scattering by 
ionized impurities. 

In the case of a bulk GaAs-like semiconductor (zinc- 
blende structure, valence band maximum at the T point), 
there are two kinds of carriers: heavy and light holes. 
Therefore, both the numerator and the denominator in 
Eq.[l3]are sums of the contributions from each of the sub- 
bands. Therefore, the resulting thermopower coefficient 
is a weighted average of the coefficients calculated sepa- 
rately for each subband at the Fermi level, which is de- 
termined by the total concentration of holes distributed 
over both subbands. Since the weights include the un- 
known parameter fiQ, which in general has a different 
value for each subband, we cannot determine this aver- 
age. However, we can consider S calculated for the heavy 
and light holes subbands separately as the limiting cases, 
as the final result has to fall within the range spanned by 
them. 

The results of computations carried out for 300 K and 
employing the 6-band k-p models are presented in Fig. 2) 
The thermoelectric power S is shown separately for the 
heavy and light holes. It is expected that the heavy holes 
dominate for the case of scattering on ionized impuri- 
ties, while the opposite is true for scattering on acoustic 
phonons (we assume that the Fermi energy is small com- 
pared to the spin-orbit splitting of the valence band) . 

The theoretical data are compared to experimental val- 
ues obtained at room temperature by Mayer et alr^ for a 
series of GaAs:Be and (Ga,Mn)As samples, in which hole 
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FIG. 4. (Color online) Room temperature thermoelectric 
power in (Ga,Mn)As and GaAs:Be as a function of hole den- 
sity changed by irradiation with high energy Ne + ions (af- 
ter Mayer et a/.— ). Lines are calculated using the standard 
six band model and GaAs parameters,— and assuming that 
either ionized impurity or acoustic phonon scattering domi- 
nates. The actual value of thermopower S should lie between 
lines obtained for heavy and light hole bands in each case. 



density was changed by irradiation with high energy Ne + 
ions that introduce compensating donor defects. We note 
that for these (Ga,Mn)As samples the content of substi- 
tutional Mn is x — 0.045, and the MIT occurs at the hole 
concentration p c w 3 x 10 20 cm" 3 ^ a value about two 
orders of magnitude higher than p c for GaAs: Be. Nev- 
ertheless, the difference between the magnitude of S in 
(Ga,Mn)As and GaAs:Be is only slight, which shows that 
the MIT has actually a little effect on the thermodynamic 
DOS, as could be expected for the Anderson-Mott local- 
ization. 

We expect that acoustic phonon scattering, for which 
r = —1/2, is relevant in GaAs:Be at 300 K. In contrast, 
in the case of (Ga,Mn)As, owing to the proximity to the 
MIT, the mobility is expected to increase with the carrier 
energy. At the same time, an additional compensation 
by interstitial Mn makes ionized impurity scattering, for 
which r = 3/2, more significant. Thus, we can conclude 
that the data for (Ga,Mn)As fall in the range expected 
for hole transport in the GaAs valence band. 



V. LOW TEMPERATURE ELECTRONIC 
SPECIFIC HEAT 

In doped semiconductors on the metal side of 
the Anderson-Mott localization, the magnitude of 
the Fermi liquid parameter is relatively small r s — 
(47rp/3) -1 / 3 ja* B < 2.4, so that Landau's renormalization 
of the specific heat and thermodynamic DOS by carrier- 



FIG. 5. (Color online) Thermodynamic density of states as a 
function of the hole concentration in (Ga,Mn)As for various 
values of the parameter Bg characterizing spin splitting of the 
hole subbands. The corresponding values of the Sommerfeld 
electronic specific heat coefficient 7 are shown on the right 
axis. 



carrier interactions is of a minor quantitative importance. 
Hence, we compute the electronic specific heat in the 
low-temperature limit according to Cy = jT, with the 
Sommerfeld constant 



I 2 

7 = y «bPf, 



(16) 



where pp = dp/dEp is determined for holes in the va- 
lence band with no carrier-carrier interactions taken into 
account. We plot in Fig.[5]the dependence of p-p and 7 on 
the hole concentration p for (Ga,Mn) As with various val- 
ues of the spin splitting parameter Bq = Apj3M / {Qgp-o). 
The numerical values for the samples A and B of Yulda- 
shev et al. are 0.37 and 0.39 mJ mol -1 K -2 , respectively, 
hence the values of 7X0 are small compared to the crit- 
ical anomaly of the specific heat, validating the results 
shown in Fig. [T] 



VI. LOW TEMPERATURE CONDUCTIVITY 

Near the MIT, temperature and magnetic field de- 
pendencies of conductivity <j(T, H) are determined 
by quantum phenomena specific to Anderson-Mott 
localization^^ These striking effects result from single- 
particle interferences of scattered waves and/or from 
scattering-driven interferences of carrier-carrier interac- 
tion amplitudes. It was suggested within this frame- 
work that the magnitude of conductivity changes in 
(Ga,Mn)As at low temperatures points to the value of 
DOS actually expected for the GaAs valence band. 26 
More recently, Neumaier et alr^ carried out compre- 
hensive studies of conductance in various dimensionality 
structures of ferromagnetic (Ga,Mn)As at subkelvin tem- 
peratures. Since the external magnetic field has no effect 
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on <j(T) in this regime^ the single-particle Anderson 
localization term, destroyed presumably by the demag- 
netizing field, does not contribute to cr{T) in this ferro- 
magnetic semiconductor below 1 K. At the same time, 
the study of <r(T) at the dimensional crossover 2D — > 3D 
allowed to determine the magnitude of the diffusion coef- 
ficient D = a I (e 2 p-p)^ In this way the value of DOS at 
the Fermi level for charge excitations pp was determined 
and found to be slightly smaller than the one expected 
for the GaAs valence band.— ^ 

To supplement the above analysis, we consider o~(T) 
within the universality class for which the transport pro- 
ceeds in two subbands, whose splitting is much larger 
than k^T but much smaller than the Fermi energy. Fur- 
thermore, we make use of the value of the Landau param- 
eter describing the correlation-induced enhancement of 
the carrier spin susceptibility and the Curie temperature 
A-p = i 2 i 15 i 32 i 33 Our goal is to describe the magnitudes 
of parameters a characteri zing the rate of change of o~(T) 
displayed in Fig. 1 of Ref. [l9| for various dimensionality 
systems. In the 2D case a = F^d^ 1 I '(nht) , where accord- 
ing to theory developed for a simple isotropic band, in the 
presence of a sizable spin splitting, F 2D = (l — F/4) In 10, 
where F = 2(Ap — 1)4>2£ The theoretical value obtained 
for the sample thickness t = 42 nm, a = 1.6 x 10 7 e 2 //im 
is seen to be in a good agreement with the experimental 
finding 1.8 x W 7 e 2 /hmM 

We now apply the same procedure to the ID and 
3D cases, d = 1 and 3^ where a ~ D x ~ d l 2 . Here, 
F 1D = 1.37 and F 3D = 1.04 are anticipated theo- 
retically for F = 0.4. 4 For the experimental values 
a = 6.95 x 10 8 e 2 //im and 3.95 x 10 8 e 2 //im as well 
as for p F = 1.98 x 10 46 1/Jm 3 and 1.45 x 10 46 1/Jm 3 , 
as expected for disorder-free (Ga,Mn)As valence band 
at the hole concentrations in question , 16 ' 19 we obtain: 
a = -3.1 x 10 7 (e 2 //i)(K 1/2 /m) for d = 1 and 0.93 x 
10 7 e 2 /(/imK 1 / 2 ) for d = 3, which are close to the ex- 
perimental values a = —2.5 x 10 7 (e 2 //i)(K 1 ^ 2 /m) and 
1.5 x 10 7 e 2 /(/imK 1 / 2 ), respectively, shown in Figs, lb 
and Id of Neumaier et al. 19 As seen, this approach sug- 
gests a slightly higher DOS comparing to that expected 
for the GaAs-like valence band. 

It is worth noting that other authors ; 41 ' 42 analyzing 
cr(T) up to 4 K, found that a(T) = ct + AT a , where 
a = 1/3. This dependence was interpreted in terms of a 
renormalization group equation 4 - applicable close to the 
MIT, where a < AT a and then 1/3 < a < 1/2 in the 3D 



casei ' Furthermore, the apparent value of a can be re- 
duced above 1 K by a cross-over to the regime, where the 
effect of scattering by magnetic excitations onto quantum 
corrections to conductivity becomes significant 2& 

VII. CONCLUSIONS 

We have considered the thermodynamic and thermo- 
electric properties of (Ga,Mn) As including the critical be- 
havior of specific heat and magnetization, high temper- 
ature thermoelectric power, low temperature electronic 
specific heat, and low temperature conductivity. The 
available experimental data 17-19 are consistent with the 
p-d Zener model in which the carriers reside in a GaAs- 
like valence band. In particular, the critical behavior of 
specific heati 7 - can be reasonably well described assum- 
ing Gaussian fluctuations of magnetization. The magni- 
tudes of the thermoelectric power at room temperature 18 
are consistent with the theoretical results for scatter- 
ing mechanisms expected to limit the magnitude of hole 
mobility. The data for low-temperature electronic spe- 
cific heat are provided in order to stimulate correspond- 
ing experimental investigations. Finally, the magnitudes 
of density of states obtained assuming a valence band 
model are supported by the temperature dependence of 
conductivity at subkelvin temperatures, determined by 
disorder-modified carrier-carrier interaction effects, sig- 
nificant near the metal- insulator transition, 19 The prox- 
imity of this transition results also in a reduction of the 
volume occupied by spins contributing to the ferromag- 
netic order. We conclude that there is no experimental 
evidence for the enhanced density of states at the Fermi 
level, expected within the impurity-band models of fer- 
romagnetism in (Ga,Mn)As. Instead, however, the prop- 
erties of this ferromagnet are strongly affected by hole 
localization effects. 
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